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We show how a two dimensional Bose-Einstein Condensate trapped in a nonequilibrium state 
can be driven into the Rayleigh- Taylor instability if an outward in-plane force is exerted on it. If 
the condensate is inside a semiconductor, the above force could arise from an inhomogeneous strain 
applied on the host semiconductor crystal. On the other hand, for a BEC of alkali atoms, the force 
could be exerted by an inhomogeneous magnetic field turned on after the BEC has been released 
from its magnetic trap. During its expansion, the condensate will break into droplets each a separate 
BEC. Therefore, one can create BEC droplets out of a single large BEC and study the coherence 
properties of the droplets with respect to each other. The discussion on the linear onset of the 
Rayleigh- Taylor instability in a BEC is generalized to three dimensions. 



i 

C 

o 
o 



> 

o 

(N 

in 
o 
-i— > 



i 

C 

O 

o 



I. INTRODUCTION 

The Rayleigh- Taylor instability ( RT ) 0,0] is observed 
in fingers that develop on industrial smoke spreading 
from factory chimneys 0, blast remnants of supernovae 
[U, superfluid droplets from charged superffuid jets 0, 
and in exciton droplets in bulk Ge . It occurs in inho- 
mogeneous fluids with density gradients when a force acts 
in the direction of less density. Layers of constant density 
become unstable towards the formation of undulations or 
wiggles. Since a wiggle displaces mass in the direction of 
the force, a constant density layer which has undulations 
has a lower potential energy than one which is straight 
0> 0| • The wavelength of the wiggles is bounded from 
below by the surface tension or the viscosity, since there 
is a greater cost in surface energy for a wiggle of shorter 
wavelength and therefore, greater surface area, while vis- 
cosity does not permit adjacent fluid elements to have 
very different velocities as is needed for short wavelength 
undulations. A necessary condition for the RT instability 
to evolve in a fluid is that it should take a shorter time for 
sound to propagate through the medium than the time 
in which the instability grows. Therefore, s > A/t must 
hold, where s is the sound speed in the fluid, while A and 
r are the wavelength and growth time of the instability. 

Ever since BECs were first created in the laboratory, 
experiments in BECs, particularly on their collapse and 
revival , the quantum interference of two BECs and 
the superfluid to Mott-Insulator transition in optical lat- 
tices [10] have justified the hope that the laws of quan- 
tum mechanics can be studied in such systems. BECs 
made out of atoms are trapped in magnetic traps and are 
physically isolated from the walls of the vacuum chamber. 
This limits decoherence processes to interactions with un- 
condensed particles which remain in the vicinity of the 
condensate. On the other hand, BECs in semiconduc- 
tor systems usually have a shorter decoherence time, r^, 
compared to their atomic counterparts because of greater 
interaction with the host, since generally ~ r e , where 
r e is the interaction time [ll[. BECs in semiconductors 
trapped in more than one dimension are truly extended 
systems with the host environment making quantum me- 



chanical measurements of the density of the condensate 
wavefunction at different parts of the condensate, hence 
quickly destroying the long-range order. The onset of 
decoherence through the destruction of long range order 
as above is related to the Kibble-Zurek mechanism [l2l |. 
which is important in cosmology. The presence of ad- 
ditional channels of decoherence within a solid due to a 
greater interaction with the host makes it more difficult 
to make BECs in solids that can survive long enough to 
be observed. Magnon condensation has been reported 
in the magnetic compounds, Cs2 Cu CI4 and Tl Cu CI3 
[l3| . A solid state system, which is an interesting but not 
yet proven as a candidate system for Bose-Einstein Con- 
densation is a cold excitonic cloud trapped in a coupled 
quantum well (CQW) [3, M, GJ] • 

When it comes to BECs, hydrodynamic phenomena 
play a very important role. For example, vortices in ro- 
tating BECs are essentially hydrodynamic in nature. To 
make safe inferences about the quantum mechanics, it is 
necessary to have an understanding of hydrodynamic in- 
stabilities, such as the RT instability, that may develop 
in them in certain experimental conditions. Further- 
more, hydrodynamic instabilities may give rise to new 
situations in which quantum mechanics can be further 
tested. In this paper, we show how in certain conditions 
a Rayleigh- Taylor instability may develop in a BEC giv- 
ing rise to the formation of droplets. These droplets, 
each of which is a BEC, can have interesting properties 
that need future study. One interesting possibility is to 
make clones of a parent BEC and then make quantum 
mechanical operations on each of them separately. 

The Rayleigh- Taylor instability in BECs has not yet 
been investigated. This paper aims to address that is- 
sue. In Sec. II, we consider a two dimensional BEC that 
can be unstable towards RT instability. Later, we shall 
see that our results can easily be generalized to three di- 
mensions. In Sec. Ill, we describe the linear onset of 
the RT instability. In Sec. IV, we remark on the future 
evolution of the instability. In Sec. V, we discuss two 
possible BEC systems in which we might see the forma- 
tion of Rayleigh- Taylor droplets. 
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II. THE RAYLEIGH-TAYLOR INSTABILITY IN 
A BEC 

A. Initial State of the 2D BEC at T=0 

We consider a two dimensional BEC of excitons at rest 
in a harmonic trap of frequency, w. The initial potential 
is V(r) ~ w 2 r 2 /2. The Thomas-Fermi (TF) density of 
the exciton cloud, p a ~ [A — r 2 ] 1 / 2 at r < ^/(h/rnw), 
with A — (2p/mw 2 ) and the chemical potential, p. 

We now consider the release of the trap and in its place, 
the immediate application of a radial force outward from 
the former trap's center. We call this as time, t = 0. 
The condensate is in non-equilibrium at t — 0. We may 
neglect the curvature of the expanding BEC as long as 
the wavelength of the instability is much smaller than 
the radius at which it sets in. To describe our system, 
we choose the y axis as being perpendicular to the sur- 
face of the condensate, with the surface at y = 0. We 
assume that the cloud expands in the y direction, and 
that the instability develops in the s-plane in the neigh- 
borhood of y = 0. We assume that the density of the 
condensate is linear at the surface where the instability 
happens. In other words, po = pi = constant at y < —5, 
Pa = — y/S) for -8 < y < 8, and p a = p 2 = con- 
stant at y > 8. We are studying the case, p\ > pi. 
We also assume that the density varies slowly over time, 
d\npo/d t -c 1/t, where r is the time of growth of the 
instability. This is true if s » A/r, in which case the 
instability evolves slower than the time it takes sound to 
propagate [l7| . The above condition involving the speed 
of evolution of the instability is necessary for an RT in- 
stability to develop in any system. 



gradient. We shall not consider the excitations in the 
system, but will only study the hydrodynamic flow of the 
part that remains a condensate. We note that in a force 
field, the lowest energy state is not the zero momentum 
state. 

We now consider the linear stability of the conden- 
sate against perturbations in which the fluid is locally 
displaced infinitesimally by Sf(x,y,t). We describe the 
disturbance by local changes in the density, p(x, y, t) = 
Po(y, t) + 8p(x, y, t), the pressure, P(x, y, t) = P {y, t) + 
SP(x,y,t), and the velocity, u(x,y,t) = 8u(x,y,t) — 
d8r(x,y,t)/dt. The BEC behaves locally as an incom- 
pressible fluid, and we may assume that the density p of a 
given volume element of the condensate remains constant 
as the element is displaced by 8f(i). 



p(f,t + St) = p Q (f-5r(t),t) 

= Po(r,t) - Wp ■ 8r(t), 

so that the perturbation of the density is given by 
The continuity equation, 



(3) 



| + v.(^) = o. 



when linearized becomes 
dSp 



Of 



+ V • {p 5u) = 0; 



(4) 



(5) 



(6) 



leading to 



Onset of the Rayleigh- Taylor instability in a 
BEC at T=0 



The hydrodynamic equation of motion for a Bose- 
Einstcin condensate at T=0 is, 



dv _ VP _ h 2 , 
dt mp 2m 2 



1 y/p J 



m V 2 



(1) 



where U, the external potential, includes both trap and 
applied potentials, P — pp is the pressure, with p, the 
chemical potential. When all the forces are perfectly bal- 
anced, the condensate is in static equilibrium, 



2m 



{^} +mpo g = 0, (2) 



VPn 



where the subscript denotes unperturbed quantities and 
mg = — VJ7 is the total external force on each particle. 
Then, at t = 0, the trap is released and immediately 
a radial outward force is applied. Therefore, from t — 
onwards, the condensate is in non-equilibrium. The 
condensate starts expanding due to the new potential 



V • 5tf = 0. 



(7) 



The linearized equation of motion for the condensate with 
only the lowest order terms in the perturbation is, 

dSu Ii> 2 
mp — = -V8P + m5pg-—\J{\J 2 8p). (8) 
at 4m 

The line of discontinuity in density, po, will become 
slightly deformed in the perturbed state. Following 0], 
let the line be then defined by y s + Sy 8 , with 



Sy s (x,t) = 8y s (0)e i{kx - wt) 



(9) 



. The discontinuity in normal stress required by equilib- 
rium is T s d 2 8y s /dx 2 , where T s is the 'effective tension' 
of the condensate, which we define as the kinetic energy 
per length along the locus of the density drop, similar to 
the surface tension defined in Ref. \}M . 



h 2 r s/2 

T, - ^- / dy\V^\ 2 



2m 



6/2 



ft 2 (P1-P2) 2 

4m<5 (pi + p 2 ) 



(10) 



assuming a linear fall in the density at the defect over 
a surface width, 8. Therefore, we should add a term 



3 



corresponding to the effective tension to the right side of 
Eq. ijHJ), which leads to, 



dSu T-,cr, c * h T-, /t-,2 <r \ 

mpQ—-—= — voP + mgdpy — - — V(V op) 
ot Am 

- yk 2 TT s Su y S{y), 



(11) 



where we have used dSy s /dt = Su y (y = 0), and therefore, 
6y s = rSu y (y = 0). 

Following Chandrashekhar 0] , we look for an instabil- 
ity of transverse wavelength k such that any perturbed 
quantity, Q, grows in time as, 



5Q(x,y,t) = Q'{yy k * +t /\ 
Taking the x-component of Eq. (jTTJ) , 

t Am 



Since, from Eq. (g]), 



dp 
dy 



Sp = —T—j—Suy, 



(12) 



(13) 



(14) 



is nonzero only at y = 0, solving Eq. (fTB")) for y ^ gives 
mpoSu x 



= -ikSP. 

T 

Similarly the y-component of Eq. (|11[) . for y ^= 0, is, 

mpoSu y dSP 
t dy ' 

Differentiating Eq. ([15]) with respect to x leads to 

mpo d5u x _ 2 
t dx 

Eqs. © and (fTTJ) lead to, 



(15) 



(16) 



(17) 



mpo dSuy ( 2 



dy 



k z 6P. 



(18) 



Differentiating with respect to y, and using Eq. (|16p . we 
find. 



mp d 2 Su y , 2 d5P mp k 2 

— k — — = ou y . (19) 



t dy 2 



Therefore, 



( — 

\dy 2 



dy 



k z \8u v = 0. 



(20) 



The condensate is unstable to velocity perturbations of 
the form 8u y = Ae 1 ! 1 e ±k v+ lkx , With boundary condi- 
tions 5u y , d5uy/dy — > as y — > ±oo, the general solution 
of Eq. (T21I)) with an interface at ?/ = is, 



fey = Ae ky , y < 



(21) 



with k > 0, since 5u y must be continuous across the 
interface. It should be noted that 8u x need not be since 
the condensate is inviscid. 

To determine the solution and thus the growth of the 
instability, we require the boundary conditions across the 
interface. Taking the y-component of Eq. (fTl"]) , we find, 



dSP 
dy 



mpoSu 



v 



T 

Th 2 d 



rmgSuy-j^- — k 2 TT s Su y 5(y) 
uy 



+ 



4m dy I V dy 
Integrating Eq. (T2"2"]) from y — 0— to 0+ gives, 



(22) 



A[6P] 



-A [mgrSiiyp] — k 2 TT s 5u y (y = 0) 



Am 



-A 



v y-dy-^ 



(23) 



where A denotes the difference across the interface. On 
the other hand, Eqs. 0, (13]) and {H} lead to, 



dy 



■po 



(24) 



since c?/?o / dy 



for y 7^ 0, and therefore, 
yuyu/uyj) — 0. Eliminating A(SP) from Eqs. 
23} and (|24|) . we find the boundary condition, 



A[V 2 {uydp /dy)} 



A 



m duy 
k 2 T dy 



■fin = A [mgrSuyiy = 0)p Q }+k 2 TT s Su y (y = 0), 

(25) 

which, after substitution for Su y using Eq. (|21|) . and Eq. 
(Unj), gives the relation between the growth time of the 
instability and the wavevector, 



1 



k l/2 



( PX ~ P2 \ 
\pl + Pi' 



h 2 k 2 /pi 

Am 2 5 V pi 



P2 



P2 



-1/2 



(26) 



For p\ > p2j the solution is unstable (r is real and pos- 

r i 1 / 2 

itive) for < k < k c where k c = < m(pi — pi)gjT s \ 

Using dr/dk — in Eq. (I26|) gives us the fastest growing 
mode, which has a wavevector and growth time, 



fc c /V3; 



3 3 / 2 h 



\Am{g 3 Sy 



C. Future Evolution of the RT Instability 

We expect the instability to progress into the nonlin- 
ear regime through the development of spikes and bub- 
bles [a, 0]- Spikes longer than their diameters are un- 
stable to varicose perturbations leading to their breakup 
into droplets. Spike breakup can also occur due to the 
Kelvin-Helmholtz Instability at the edges of the spikes 
in both 2D and 3D cases. Spike breakup leads to droplet 
formation in both cases instead of a mixing zone when 
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surface tension is significant, since droplets minimize sur- 
face energy. The diameter of the droplets is then half the 
wavelength of the fastest growing mode, A/2. 

For low non-zero temperatures, the motion of the con- 
densate and the excitations are decoupled. Since the ex- 
citations are dilute and of long wavelength, the exciton 
fluid continues to be RT unstable with the droplet size 
determined by the effective tension of the condensate and 
the applied stress or electric field gradient. For weak in- 
teractions and low temperatures, PoUq << T « T c , the 
condensate moves inviscidly in a mean field due to the 
excitations in the Hartree-Fock approximation. If the ex- 
ternal potential is modified to include the mean field due 
to excitations as a correction, the discussion in Sec. Ill 
can be used to show that the BEC is RT unstable with 
the applied force and the effective tension again deter- 
mining the droplet size. At T > T c , however, when the 
fluid is normal and has viscosity, the character of the flow 
changes. Though the fluid can still be RT unstable, the 
droplet size is now determined by the viscosity and ap- 
plied force. Therefore, one can observe in experiments 
the variation in droplet size as T is increased beyond T c . 

III. APPLICATION TO SPECIAL CASES 

Below we investigate the effects of the Rayleigh- Taylor 
instability in two systems of special interest. 

A. Atomic BECs 



tion of g ~ 4.55 — 455 cm/s 2 , at a distance of 2 pm away 
from the trap center. A BEC contained in a trap of fre- 
quency, 24 Hz will have a radius ~ 2/im. After the BEC 
is released from this trap, an inhomogeneous magnetic 
field can be turned on so that the force exerted points 
outwards from the center instead of inwards. The BEC, 
which is now in a nonequilibrium state will then begin 
to expand. If the force of expansion is of the appropri- 
ate magnitude, there will be an RT instability driven by 
the force of expansion which will lead to the formation 
of RT droplets. A strong outward force corresponding 
to an acceleration, g ~ 5 x 10 5 cm/s 2 , leads to a fastest 
growing wavelength of 0.1 pm (Eq. |27|) . if we assume 
that the surface width, 6, is almost the same in magni- 
tude as the wavelength. This leads to droplets of size, 
> 0.05 fim, which can be observed. The growth time 
for such a droplet is given by Eq. [27] to be f ~ 3 /is. 
Weaker forces leads to larger wavelengths for the insta- 
bility. Since for an observable instability, we need the 
wavelength, A < irD, where D is the diameter of the 
cloud, we find that the maximum value of acceleration at 
the surface is, g m ax = 0.8 cm/s 2 , which leads to a mini- 
mum growth time for observable instability, f m i n ~ 0.02 
sec for X m ax — D ~ 4 pm (Eq. |27) . 

A suitably designed experiment can therefore be used 
to create Rayleigh- Taylor droplets in an atomic BEC. In 
an actual experiment, the condensate will be allowed to 
expand to several times its initial size before measure- 
ments are made. Therefore, the droplets seen at mea- 
surement will be scaled up. 



Consider a typical experiment in which we have a BEC 
of alkali atoms in a three dimensional magnetic trap. 
When the trap potential is switched off, the condensate is 
in a non-equilibrium state. Repulsive interactions make 
the condensate expand slowly. After a a time of flight of 
several nanoseconds, a snapshot of the freely expanding 
BEC is usually taken. We propose one modification of the 
experiment which can lead to the formation of droplets of 
BEC at the surface of the condensate while the whole gas 
expands. We suggest that after turning off the magnetic 
trap, we turn on an "inverted magnetic trap" with a po- 
tential profile, V(r) ~ —1/2 w 2 r 2 , resulting in a radial 
outward force acting on the expanding condensate. This 
outward force will result in a Rayleigh- Taylor instability 
at the surface of the expanding condensate. 

It is simple to generalize to three dimensions the results 
of Sec. III. We continue to choose y as the linear axis 
perpendicular to the surface. In 3D, the Eqs. (TTT]) . (f2"D|) - 
and the later expressions for the growth time and 
wavelength, Eqs. (j26|) and (|27| hold exactly with k now 
being given by fc 2 = fc 2 + fc 2 + fc 2 , and every perturbed 
quantity, Q, growing in time as, 



SQ(x,y,t) = Q'( y )e i *- iB + tt »»+ ifc '*+*A 



(28) 



B. Excitonic BEC in Layered Semiconductors? 

The possibility exists that an indirect exciton cloud in 
a layered semiconductor if cooled below Tc may form two 
dimensional BECs trapped in impurity driven in-plane 
shallow traps [3, Oil- This type of BEC has not yet been 
conclusively observed. If, however, such a BEC can be 
created in the future, it may be possible to overturn the 
trap by applying an outward radial force through an inho- 
mogeneous stress. This would lead to a two-dimensional 
BEC in a non-equilibrium state, which is driven to ex- 
pand by the radial force resulting in an RT instability. 
We note here that inhomogeneous stress has previously 
been used to create a harmonic trap in a semiconductor 
bilayer [2(|. It is reasonable to assume that it may be 
possible to use the same means to create an outward, 
instead of inward, radial force field. 

For indirect excitons, m = 0.1m e , p ~ 10 16 c™ ' 
the wells, and ms 2 ~ pUo ~ 0T — 1 meV 14|. If g 



For 87 Rb atoms, a magnetic trap of frequency, 24 — 240 
Hz can exert a restoring force equivalent to an accelera- 



cm m 
10 15 

cm s~ 2 as achieved in Ref. [2fJ through inhomogeneous 
stress and electric fields, and making the reasonable as- 
sumption that S ~ A, under which conditions the RT 
instability still exists, we find A ~ 5.4 pm and f ~ 0.4 
ns. Moreover, it can be verified that s > X/f. Therefore, 
we would expect BEC droplets, 5-10 microns in diame- 
ter, but as the exciton cloud starts to expand radially, 
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the size of the droplets will scale accordingly. 

IV. CONCLUSIONS 

While different statistical distributions are used to cor- 
relate data on the drop size in normal fluids, such as 
the Nukiyama-Tanasawa law 0, l2ll ] , which is arrived at 
through assumptions on droplet formation, a theory for 
the distribution in the droplet size in a quantum fluid 



such as a BEC needs to be developed. We also need 
new experiments to study BEC droplets to further our 
understanding of the quantum mechanics. 
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